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SECTION  1.  GENERAL.  DESCRIPTION 


1.1  INTRODUCTION 

Present  methods  for  predicting  internal  wave  fields  produced  by 
objects  moving  in  stratified  media  are  limited  in  their  physical  realism, 
flexibility,  and  economy.  Recent  efforts  to  define  standards  for  internal- 
wave  prediction  have  helped  to  clarify  the  shortcomings  of  existing 
techniques,  which  occur  in  one  or  more  of  the  following  areas: 

©  Modelling  the  dynamics  of  turbulent  wake  *;owth,  buoyancy  trans¬ 
port,  and  wake  collapse. 

o  Modelling  the  propagation  of  radiated  internal  waves  in  a  medium 
of  arbitrary  buoyancy  stratification. 

•  Maintaining  numerical  precision  and  stability  in  a  field  calcula¬ 
tion  encompassing  a  wide  range  of  horizontal  and  vertical  length 
scales,  at  acceptable  computing  cost  and  speed. 

Some  of  these  shortcomings  are  fundamental  and  unavoidable,  while 
some  are  accidental,  depending  on  the  approach  used.  For  example,  a 
three-dimensional  finite-difference  calculation,  superior  for  near-field 
flow  and  wake  dynamics,  is  severely  limited  in  the  upper  length  and  time 
scales  it  can  achieve,  simply  because  of  limits  on  available  computer 
storage,  speed,  and  accuracy.  On  the  other  hand,  limitations  in  existing 
analytic  approaches  can  be  overcome.  One  such  approach  now  in  use  by 
TRW  employs  a  three-layer  model  of  density  stratification  that  is  too 
simplified  to  reproduce  wave  propagation  on  real  thermoc.lines ,  although 
it  is  perfectly  adequate  for  its  original  purpose,  which  was  to  provide 
estimates  of  field  magnitudes  and  shapes  [2], 
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A  more  general  analytic  cede  can  be  envisioned,  one  based  on  a 
Fourier /normal-mode  expansion  of  the  linearized  field  equations  for  an 
arbitrary  profile  of  stratified  buoyancy,  that  would  meet  the  require¬ 
ments  for  generality,  scale  and  efficiency.  Like  its  predecessors, 
however,  such  a  code  would  have  to  depend  on  parametric  source  models  to 
represent  the  excitation  of  internal  waves  by  the  moving  body  and  col¬ 
lapsing  wake,  source  models  that  will  have  to  be  devised  and  validated 
elsewhere . 

A  two-fold  approach  has  been  suggested  [3]  in  which  ah  efficient  code 
of  the  analytic  type  is  developed  for  routine  simulations,  while  at  the 
same  time  a  "research"  code  of  the  finite-difference  type  is  developed 
to  provide  validated  source  models  for  the  analytic  code. 

This  two-volume  report  describes  a;  prototype  analytic  code,  XMODE, 
which  has  been  developed  at  PDA  to  provide  inexpensive  simulations  of 
radiated  internal-wave  fields  in  a  variety  of  realistic  thermoclines. 

XMODE  represents  an  improvement  over  existing  codes  in  two  important 
respects:  it  contains  an  efficient  eigenfunction  generator  for  the  normal¬ 
mode  analysis  of  an  arbitrary  input  density  prpfile,  and  it  abandons  the 
stationary-phase  method  of  field  calculation  in  favor  of  a  more  nearly 
exact,  and  very  rapid,  inversion  in  rectangular  coordinates  via  the  Fast 
Fourier  Transform  [4j. 

Volume  I  contains  a  general  description  of  XMODE  and  a  detailed 
derivation  of  the  algorithm  for  representative  body  and  wake  sources. 
Volume  II  is  a  user's  guide  to  XMODE,  containing  a  description  of  the 
computer  routines,  operating  instructions,  and  FORTRAN  listings. 


Some  of  the  procedures  and  computer  subroutines,  those  having  to 
do  with  eigenfunction  generation,  are  identical  to  those  in  the  computer 
code  ZMODE,  which  is  an  integral  part  of  XMODE  and  which  has  been  described 
fully  in  an  earlier  report  [1].  This  material  is  not  duplicated  hare; 
potential  users  Gf  XMODE  should  obtain  Reference  1  as  a  companion 
volume  to  the  present  report. 

1.2  OUTLINE  OF  METHOD 

The  algorithm  used  by  XMODE  is  based  on  a  Fourier/normal-mode 
expansion  of  che  linearized  equations  of  stratified  flow,  with  simple 
local  sources.  As  in  other  formulations  the  transformed  solution  is 
obtained  as  an  algebraic  combination  of  eigenfunction  amplitudes  and 
dispersion  quantities.  However,  the  present  formulation  is  distinct 
in  that  the  transform  is  inverted  in  rectangular  coordinates,  first  along 
the  x-direction  (track)  by  analytic  means,  then  along  the  y-direction 
(crosstrack)  by  a  numerical  Fast  Fourier  Transform  (FFT).  The  partial 
crosstrack  transforms  obtained  in  the  first  step  for  each  normal  mode 
are  complex  functions  whose  amplitudes  are  invariant  with  respect  to 
x  and  whose  complex  phases  are  linearly  proportional  to  x.  This  property 
allows  the  partial  transforms  to  oe  assembled  rapidly  from  independent 
amplitude  and  phase  factors  computed  and  stored  ahead  of  time.  At  each 
value  of  x  the  numerical  phasa-and-sum  operation  is,  for  twenty  modes, 
no  more  time-consuming  than  the  FFT  computation  of  the  crosstrack  field, 
so  that  the  full  efficiency  of  the  FFT  is  realized. 

1.3  ADVANTAGES  AND  LIMITATIONS 

The  logic  of  XMODE  is  designed  for  efficient  computation  of  radiated 
internal  wave  field  properties  at  the  ocean  surface,  or  on  a  horizontal 
plane  (x-y)  at  any  specified  depth.  Computation  of  the  field  on  the  ver¬ 
tical  crossplane  (y-z)  is  somewhat  less  efficient,  and  confutation  on  a 
vertical  plane  (x-z)  is  least  efficient  of  all.  All  points  in  the  plane 
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are  computed  with  equal  precision;  unlike  previously  used  stationary- 
phase  methods,  the  FFT  in  rectangular  coordinates  has  no  difficulty  with 
points  on  or  near  the  coordinate  ares  (x  =  0  and  y  =  0).  originating  on 
the  source. 

Other  specific  assumptions  and  features  are: 

«  Linearized  fluid  equations  with  dipole  body  and  quadrupole  wake 
sources. 

0  Choice  of  field  quantity  among  vertical  displacement,  scalar 
strain;  x-  and  y-  components  of  velocity,  and  scalar  strain 
rate. 

®  No  ambient  shear  flow. 

£  I-lo  vertical  momentum  in  wake. 

&  Source  speed  must  exceed  maximum  ambient  phase  speed  (Super- 
Froude  source) . 


SECTION  2.  SPECIFICATIONS  OF  THE  XMODE  CODE  AND 
SAMPLE  CALCULATIONS 

2.1  RESOLUTION  AND  DYNAMIC  RANGE 

Each  crosstrack  field  contains  256  physical  resolution  elements  Ay, 
field  half-width  =  256  •  Ay; 

the  width  and  resolution  size  are  somewhat  adjustable,  but  Ay,  which  should 
approximate  the  reciprocal  of  the  upper  wavenumber  limit  in  the  ZMODE  cal¬ 
culation,  is  constrained  by  numerical  stability  considerations  to 

Ay>  Za/15  , 

where  Z  is  the  total  depth  of  the  active  thermocline.  For  Z  =  300  m 

3  cl 

this  permits  a  crosstrack  resolution  of  20  m  and  a  field  half-width  of 
5120  m.  The  number  of  points  calculated  per  crosstrack  field  is  optional 
among  256,  512,  1024,  and  2048,  for  redundant  sampling  at  multiples  of 
2,  4,  and  8  times  the  Nyquist  frequency  Ay 

Any  number  of  equally-spaced  crosstrack  fields  can  be  requested  in 
a  given  calculation,  with  an  arbitrary  track  spacing  Ax. 

The  maximum  number  of  normal  modes  is  20,  so  that  the  effective 
vertical  resolution  is  AZ  -  2/20. 

C7. 

Since  the  source  transforms  are  generated  directly  from  analytic 

formulas  and  stored  in  floating-point  format  the  effective  dynamic  range 

of  the  tabulated  transforms  is  very  large,  limited  only  by  the  eigenmode 

-20 

convergence  precision  of  ZMODE,  which  is  in  the  neighborhood  of  10  , 

or  double  the  eigenvalue  c r  .vergence  precision. 
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Errors  introduced  by  the  finite,  discrete  FFT  are  of  two  simple  types:, 

accumulated  roundoff,  which  is  a  few  times  machine  precision,  and  aliasing, 

which  superposes  fields  from  image  sources  spaced  at  multiples  of  the  field 

width  from  che  actual  track.  The  image  fields  are  negligible  so  long  as 

the  simulation  time  is  short  cc  tared  to  the  time  required  for  the  fastest 

propagating  field  components  to  traverse  a  half-width;  when  this  artificial 

periodicity  is  taken  into  account,  the  images  will  not  affect  estimates  of  the 

-13 

spectral  content,  which  can  be  recovered  with  a  precision  of  10 

2.2  CORE  REQUIREMENTS  AND  EXECUTION  TIMES 

T1  XMODE  source  deck  contains  1120  cards.  On  the  CDC  7600,  both 
the  RUN  compiler  and  the  FTN  optimizing  compiler  produce  an  object  code 
occupying  less  than  60,000  words  (162,000  octal). 


The  execution  tirors  are  proportional  to  the  number  of  modes  requested 
in  a  calculation.  The  values  listed  below  are  for  the  full  20  modes. 

Execution  Times,  Seconds 


Procedure 


•  Dispersion  Tables 
and  Mode  Amplitudes 

»  Recalculate  with 
new  Source  Depth 

•  Crosstrack  Source 
Transforms 

®  With  Printer  Plots 

a  Radiated  field  on 
100 (x)  by  1024 (y) 
array 


Module 

RUN 

FTN 

ZMOBE 

11.2 

6.6 

ZM0DE 

4,1 

2.5 

SOURCE 

0.02 

0.01 

SOURCE 

0.30 

0.38 

XFIELD 

6.4 

5.3 

6 


'  Particularly  noteworthy  is  the  efficiency  achieved  by  the  XF1ELD 

module.  Of  the  64  milliseconds  computation  time  per  1024-point  crosstrack 
field  (RUN  compiler),  half  is  devoted  to  the  20-mode  phase-and-sum,  and 
half  to  the  FFE. 

The  computation  time  of  the  dispersion  tables,  the  most  time-consuming 
operation,  will  depend  also  on  the  chosen  number  of  eiganmode  integration 

•  I 

steps,  the  convergence  precision,  and  tabular  wavenumber  density;  the  values 
}  200,  10  and  41  used  in  the  test  case  are  felt  to  be  typical  for  deep- 

, 4  ocean  calculations.  For  a  series  of  field  calculations  with  differing  source 

speed  or  field  type,  a  single  dispersion  calculation  suffices.  A  change  in 
,  |  source  or  field  depth  requires  a  new  calculation  of  „ne  eigenfunctions  so 

that  the  mode  amplitudes  at  the  new  depths  can  be  tabulated.  When  converged 
eigenvalues  are  already  present  the  eigenvalue  search  can  be  bypassed, 
shortening  the  (RUN)  time  from  11.2  to  4.1  seconds,  as  indicated. 

;  t 

•  «  2.3  SAMPLE  CALCULATIONS 

The  sample  calculations  illustrated  below  simulate  a  source  60  meters 

•  ! 

deep  moving  at  2.0  kt  and  16  kt.  The  depth  chosen  places  the  source 
:  i  15  meters  below  the  sharply  eroded  boundary  of  the  test  thermocline, 

I  1  shown  in  Figure  2-1,  which  has  been  taken  from  data  obtained  on  the 

,  ^  oceanographic  research  vessel  FLIP.  The  corresponding  profile  of  Vaisala 

,  .  frequency  and  the  ZMODE-derived  internal  wave  dispersion  plot  are  shown 

, .  in  Figures  2-2  and  2-3. 

I I 

The  crosstrack  source  transforms  computed  by  the  SOURCE  module  are 

• 

|  available  in  numerical  and  graphical  output,  figure  2-4  is  a  sample  of  the 

numerical  output,  which  includes  transform  amplitude  versus  both  the  scalar 
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wavenumber,  k,  and  its  crosstrack  component,  k^.  The  graphical  output. 
Figures  2-5  and  2-6,  plots  the  squared  transform  amplitude  against  k^. 

These  two  examples  show  the  modal  energy  spectra  of  crosstrack  surface 
current  tor  the  first  six  modes  at  2.0  kt  and  16.0  kt  respectively.  For 
greater  clarity,  the  usual  XMODE  printer-plot  subroutine  has  been  replaced 
by  an  equivalent  pen-plot  routine  specific  to  the  RDA  plotting  hardware. 

The  low-speed  and  high-speed  cases  are  dominated  respectively  by  body 
and  wake  excitation;  note  the  flat  spectra  characteristic  of  the  latter 
and  peaked  spectra  of  the  former,  especially  in  the  first  mode,  where  the 
low-wavenumber  components,  whose  phase  speed  is  about  1.6  kt,  are  nearly 
resonant  at  2.0  kt. 

The  output  of  the  field  calculation  in  XFIELD  is  placed  directly  on 
magnetic  tape.  At  the  sa_e  time  XFIELD  generates  a  summary  printout,  as 
in  Figure  2-7,  containing  the  peak  amplitude  and  corresponding  location 
for  each  crossfield.  The  boundary  amplitudes  are  useful  at  short  simu¬ 
lation  times  as  an  indication  of  the  transform  roundoff  and  aliasing  errors. 

Figures  2-8  and  2-9  are  raster  plots  made  directly  from  the  sample 
calculation  output  tape,  depicting  crosstrack  surface  current  associated 
with  the  radiated  internal  waves  in  modes  1-20.  Plot  ranges  of  0-2.0  km  and 
0-750  m  crosstrack  are  used  for  the  2.0  kt  and  16.0  kt  cases,  respectively, 
or  40%  and  15%  of  the  calculated  field,  for  best  rendering  of  the  features. 
The  downtrack.  coordinate  is  in  time  units,  for  both  cases  spanning 
5000  sec  in  increments  of  50  sec. 
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***  SQUARED  CROSSTRACK  TRANSFORM  *** 


MOOE  3  QUANTITY  •  •  UY 


K  CY/XM 

XY  CY/KK 

BODY  CQ*m)**2 

wAKE  (Q*M)**2 

0,000 

0,000 

C, 

1,4060*06 

,500 

,988 

2,9006*05 

1,2940*06 

1,000 

,980 

7,0000*05- 

1 , 1566-06 

1.500 

1,973 

9,823E*0S 

1 ,0880-06 

2,000 

1.969 

1 , 1586-04 

1 ,0600*06 

2,500 

2,466 

1,2356*04 

1,0170-06 

3,000 

2,963 

i ,2I2E*04 

9,3150-07 

1,500 

3,461 

1,1086*04 

8,1530-07 

4,000 

3,959 

9e5?8E*05 

6,8950*07 

a, 500 

4,457 

7.939E-05 

5,6906-07 

5,000 

4,956 

6,369E®05 

4,61 6E-07 

5,500 

5,455 

4,9906-05 

3,6966-07 

6,000 

5,955 

3,B16E*05 

2,9290-07 

6,500 

6,454 

2,8786-05 

2,J02fc*07 

7,000 

6,954 

2,14HE»05 

1,7966*07 

7,500 

7,454 

1 ,577E*05 

1,3926-07 

8,000 

7,954 

1,150£*05 

1,0716*0? 

8,500 

8,455 

8,1206*06 

8,1816*08 

9,000 

8.9S5 

5,96SE*06 

6,1906-08 

9,500 

9,455 

4 ,230E*Oto 

4 ,651 E»08 

10,000 

9,956 

2,98^0-06 

3,4S2E-08 

10,500 

10,456 

2 , Q  756-06 

2,5316-08 

1 1,000 

10,957 

l  ,425E»06 

1,0296-08 

11,500 

11,457 

9,6586-07 

1,3010*09 

12,000 

11,957 

6,4410-07 

9,1006-09 

12,500 

12,456 

4,2190-07 

6,2416-09 

13,000 

12,958 

2,7086-0? 

4,1806-09 

13,500 

13,459 

1  ,6980-07 

2,7426*09 

19,000 

13,959 

l.OlbE-07 

1 , 7446*09 

19,500 

14,460 

6, 1048-08 

1,0716-09 

15,000 

14,960 

3*4440*08 

6,2960-10 

15,500 

15,461 

1  , 0330*08 

J.467E-10 

16,000 

15,961 

0,9600*09 

1,7?3E*10 

16,500 

16,462 

3,8241*09 

7,064E«11 

17,000 

16,962 

1,2660*09 

2,?036«ii 

17,500 

17,463 

2, 1606*10 

4,7926-12 

18,000 

17,963 

3.175E-12 

7,3070*14 

18,500 

18,465 

2,1636*10 

5,1656*12 

19,000 

18,964 

6,1400*10 

1,5190*11 

19,500 

19,464 

l  ,062E«09 

2,7216*11 

20,000 

19,965 

|,492E®09 

5,95TE®11 

Figure  2-4.  Sample  SOURCE  Module  Printout  for  Mode  3,  Source  Speed  2.0  Kt 
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WAVENUMBER  ky(CYCLES  KM  } 


Figure  2-6.  Excitation  Spectra  of  Surface  Crosstrack  Velocity  (Squared  Sourcp 
Transforms)  for  Modes  1-6.  Source  Speed  16.0  Kt 


*  A  * 

FIELD  OUTPUT  SUMMARY  *** 

4 

QUANTITY 

uv 

4 

SPEED,  K? 

16. 

00 

4 

DEPTH,  m 

60.0 

4 

XMA5UHID7H 

205 

TIME,  SEC 

TRACK,  M 

peak  signal 

AT  Y,  M 

BOJNDARY 

50,0 

41  i  ,5 

1.1372-04 

30,0 

1,8062-10 

ioo.o 

623.1 

1.292E-04 

30,0 

1.331E-10 

150,0 

1234,6 

1 ,196E«04 

35,0 

1 .21  IE-10 

200,0 

1646,2 

1.022E-04 

45,0 

1 , 2892-10 

250,0 

2057,7 

9.696E-05 

60,0 

1.60  IE- 10 

300,0 

2469,3 

1.144E-04 

90,0 

2,1 162-1 0 

350,6 

2660,6 

1.337E-04 

95,0 

2.787E-10 

000,0 

3292,4 

1 .48BE-04 

95,0 

3.545E-10 

050.0 

3703,9 

1.564E-04 

95,0 

4.307E-10 

500,0 

4115.4 

1 i580E-04 

100,0 

4,9962*10 

550,0 

4527.0 

1.527E-04 

100,0 

5 , 502E»1 0 

600,0 

4938,5 

1.425E-04 

105.0 

5 , 788E-1 0 

650.6 

5350,1 

1.293E-04 

105.0 

5,7922*10 

700,0 

5761,6 

1.1U3E-04 

110,0 

5.474E-10 

750 , 0 

6173,2 

9,915E»0S 

110,0 

4.794E-10 

600,0 

6584,7 

8.449E-05 

115,0 

3.710F.-10 

650,0 

6996,2 

7.747E-05 

160,0 

2, 1 74E*  1  0 

900,0 

7407,8 

7.267E-05 

160,0 

1.331E-11 

<350.0 

7819.3 

6.906E-05 

220,0 

2.4&2E-10 

1600,0 

8230.9 

6.571E-05 

225.0 

5.649E-10 

1050.0 

8642,4 

8.376E-05 

30.0 

9.453E*10 

1100,0 

9054,0 

1 .129E-04 

30,0 

1.388E-09 

1150,0 

9465,5 

1 .444E-04 

35,0 

1.894E-09 

1260,0 

9877.1 

1,775E»04 

40,0 

2 , 463E“09 

1250,0 

10288,6 

3.099E-Q4 

45,0 

3.098E-C9 

1300,0 

10700,1 

2 , 368E"04 

50,0 

3.602E-09 

1330,0 

11111.7 

2.553E-04 

50,0 

4,581E*09 

laoo.o 

11523.2 

2.683E-0U 

55,0 

5 , 44SE“09 

i«so,o 

11934,8 

2.712E-04 

60.0 

6.394E-09 

1500,0 

12346.3 

2 , 656E»04 

60.0 

7,4422*09 

1550,0 

12757,9 

2.588E-04 

65,0 

B.592E-09 

1600,0 

13169,4 

2,440E«O4 

70,0 

9.852E-09 

1650,0 

13561,0 

2 , 31 OE«06 

70,0 

S,123E"09 

1700,0 

13992,5 

2.159E-04 

70,0 

1.273E-0S 

1750,0 

14404,0 

2.043E-04 

70,0 

1.4386*03 

1800,6 

14015,6 

1.991E-04 

75,0 

1,6192*08 

1830,0 

13227,1 

2.008E-04 

75,0 

1,81 82*08 

1<900,0 

15630,7 

2.066E-04 

75,0 

2,0O0S«08 

1950,0 

16030,2 

2 , 1  44E"04 

80,0 

2,28hE-08 

2000,0 

16461,9 

2.248E-04 

80,0 

2.562E-03 

2030,0 

1687S.J 

2.385E-04 

05,0 

2«B71E*08 

2100,0 

17284,8 

2.548E-04 

90,0 

3 , 2 182-08 

2150,0 

17696,4 

2,7302*04 

95,0 

3,6102*08 

2200,0 

18107,9 

2.908E-C4 

100,0 

4 ,0332*06 

2230,0 

18519,5 

3.058E-04 

100,0 

4 ,3552*08 

2300,0 

18931,0 

3,21 2E-0U 

105,0 

5,1252*08 

2330,0 

19342,6 

3 ,282E°04 

110,0 

5.771 £-06 

2000,0 

19734,1 

3,3412*04 

110,0 

6.505E-08 

2450,0 

20165,7 

3.33SE-04 

115,0 

7,3352*08 

2500,0 

20577,2 

3 ,3772*04 

115,0 

3,2762*08 

Figure  2-7.  Sample  Summary  Printout  From  XFIELD  Module 


15 


P.-2074(U) 


Figure  2-8.  Raster  Plot  of  Crosstrack  Surface  Current  U  for  2.0  Kt  Source  Speed. 
Interline  Space  Represents  Amplitude  of  0.0  r  cm  sec"-. 

Track  Length  Depicted  is  6.14  Km. 


SECTION  3,  BODY-GENERATED  WAVES 


Internal  waves  are  generated  in  an  incompressible,  .^ratified  fluid 
when  a  moving  object  displaces  the  level  surfaces  from  equilibrium.  If 
the  radiated  displacement  amplitudes  are  small,  and  if  the  traverse  time 
of  the  object  is  short  compared  to  the  time  scale  associated  with  strati¬ 
fication,  then  the  excited  field  is  well  described  by  linearized  equations 
with  a  simple  set  of  singular  displacement  sources.  The  one  approximation 
that  may  be  unrealistic  for  certain  cases  is  the  neglect  of  ambient  cur¬ 
rent  shear. 

The  following  normal-mode  analysis  of  body-generated  radiation  is 
similar  to  the  general  treatment  of  Miles  [5]  .  Certain  differences  of  format 
are  dictated  by  the  rectangular  coordinates  chosen  for  ease  of  computation, 
and  these  will  be  pointed  out  an  they  are  introduced. 


EQUATIONS  OF  MOTION 


The  coordinates  me  x  (track),  y  (cross-track),  and  z  (depth,  positive 
upward,  zero  at  surface).  The  fluid  displacement  from  equilibrium  will  be 
described  by  a  vertical  compou'  at  £  and  a  twc-component  vector  5  in  the 
horizontal  plane.  In  the  usual  Boussinesq  approximation  with  small  dis¬ 
placements  and  small  departures  p  from  equilibrium  pressure,  the  equations 
of  motion  are 


£  ■*’  jHr  +  N  (z)  £ 
0  x 


(3-1) 


i  +  V'P  -  o  (3-2) 

where  V'  stands  for  'he  horizontal  gradient  (3/3x,  3/3y)  and  N  is  the 
Vaisala  frequency  associated  with  the  buoyancy  stratification 


Conservation  of  volume  will  be  modified  to  include  a  local  set  of  sources 
and  sinks  describing  the  moving  body, 

~+  VIE,  =  s(x-vt,  y,  z)  .  (3-4) 


To  eliminate  the  pressure,  we  apply  8/3z  to  (3-2),  V'  to  (3—1} ,  and  subtract, 


V'\  + 


N Vs  - 


3£ 

9z 


=  0 


(3-5) 


we  can  then  eliminate  the  §  term  by  taking  the  horizontal  divergence  of 
(3-5)  and  adding  the  result  of  (3/3z)  (3"/3t^)  on  (3-4), 


+  N2V'2? 


(3-6) 


3.2  NORMAL  MODE  SOLUTION  IN  FOURIER- .KANSFORK  COORDINATES 

Using  the  nota  ;ion  (“)  for  the  2-dimensional  Fourier  transforms  in 
the  horizontal  plane 

X  (k,z)  =  ^  e  ik’X  C  dxdy  (3-7) 

and  so  on,  with  k*x  =  k  x  +  k  y  ,  and  observing  that  every  component  of 

x  y 

a  steady  solution  emanating  from  the  moving  source  (3-4)  must  have  the 
particular  time  dependence 

r(k,z)  e  Ut  ,  u  =  kxv,  (3-8) 


we  get  from  (3-6) 


2  2  2  — 

where  k  «  k  4k  and  s  is  the  transform  of  s  at  t  °  0, 
x 


(3-9) 
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At  each  value  of  k,  the  solution  £  to  the  inhomogeneous  equation 
above  can  be  obtained  as  a  linear  combination  of  the  eigensolutions  <{> 
of  the  corresponding  homogeneous  equation 


32<j>  2  /n2 

TT  +  k  “-1!  *m  =  °’ra=0’  x>2> 
32  \  0) 

\  m 


(3-10) 


which  define  the  freely  propagating  normal  modes.  The  functions  <{> 

•k  ^ 

vanish  at  z  =  0  and  at  z  =  b,  the  ocean  bottom.  They  form  a  complete 

orthogonal  set, 


[°  2 

J  *nN  ♦»  dz  *  5nm- 

-b 


(3-11) 


over  the  vertical  interval  in  which  N(z)  is  nonvanishing.  In  these 

defining  equations  che  scalar  wavenumber  is  a  continuous  parameter  upon 

which  the  functions  <t>  and  the  eigenvalues  w  depend,  a  convention,  opposite 

m  m 

to  that  of  Miles'  treatment,  in  which  k  is  the  eigenvalue  and  w  the 
parameter.  The  present  choice  simplifies  the  logic  of  the  numerical 
algorithm. 


We  now  proceed  in  the  usual  way  to  assemble  the  eigenfunctions' into 

a  Green's  function  G,  (z,z  )  for  (3-9), 

k,o>  o  ’ 


GkiU  .  Hz  -  z0) ,  (3-12) 


Coupled  displacements  of  the  free  surface  are  very  small  and  can  be 
safely  neglected.  See  Phillips  [6]  . 
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which  in  combination  with  the  source  transform  3s/3z  will  yield  the 
displacement  transform  The  identity 


V  <j>n(z)  6  (z  )  N2(t)  =  <5(z  -  z  )  (3-13) 

/ ,  m  mo  o 

m=l 


and  the  property,  via  equation  (3-10), 


k2N2 


m 


together  imply 


(3-14) 


To  represent  the  effect  of  a  cylindrically  symmetric  body  we 
distribute  displacement  sources  along  the  submerged  track  z  =  zq,  y  =  0, 

s(x,y,z)  =  A(x)  5(y)  6 (z  -  zq) ; 


the  particular  choice 


A(x) 


const. ,  |x|  <_  L/2 
0,  [ x|  >  L/2 
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is  equivalent  to  the  Rankine  ovoid  of  volume  V  =  AL  and  dipole  length  L 
(except  for  the  effects  of  local  stratification  and  the  rigid  surface 
in  distorting  the  boundary  streamlines).  For  this  choice,  the  source 
transform  term  in  (3-9)  is 


sin  k  L/2 


3lKM>  -  v 


6' (z  -  z  ) 
o 


(3-15) 


where  the  prime  (')  is  used  here  and  subsequently  to  denote  differentia¬ 
tion  by  z.  For  cases  of  practical  interest  the  values  of  k^  =  u)m/v 
present  in  the  radiated  field  will  be  such  that 


k  L 
x 

2 


<<  IT 


so  that  the  form  factor  will  be  very  nearly  unity;  note  that  this  inequal¬ 
ity  must  hold  also  for  the  distortions  in  the  Rankine  ovoid  boundary 
streamlines  to  remain  small.  The  body  is  therefore  effectively  a  point 
monopole  for  displacement  (or  point  dipole  for  velocity) .  The  equation 
for  the  displacement  transform  then  becomes 


?  =  V  6’(z  -  zQ) 


(3-16) 


so  that  C  is  minus  V  times  the  derivative  with  respect  to  zq  of  the 
Green's  function  (3-14), 
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(3-17) 


The  subscripted  terms  in  the  above,  and  depend  on  the  scalar 

magnitude  of  k,  while  to  is  identical  to  the  track  component,  k  except 

for  the  scale  factor  v,  w  =  k  v. 

’  x 

3.3  CROSSPLANE  REPRESENTATION 


The  inversion  of  the  field  transform  £  will  be  carried  out  in  two 
steps,  first  analytically  along  the  kr  coordinate  to  produce  the  "partial 
transforms"  l, 


?(x,ky,z)  •*  2V 


1  fm  ,  -ik  'c  ,, 

J  £(k,z)  e  x  dkx, 


(3-18) 


then  along  the  k^  axis 


C(x,y,z)  =  —  j  \ (x,ky,z)  e  ikyy  dk  .  (3-19) 


The  second  step  is  carried  out  numerically  with  the  Fast  Fourier 
Transform  for  a  succession  of  x-values  at  a  given  z.  As  will  be  seen, 

'V/ 

the  analytic  form  of  the  partial  transform  £  allows  it  to  be  assembled 
rapidly  at  each  x  from  quantities  computed  and  stored  ahead  of  time. 

Each  term  in  the  modal  sum  (3-17)  for  the  field  transform  contains 
an  integrable  singularity  along  the  line  of  resonance  in  the  k  plane 
defined  by 
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(3-20) 


k  v  =  a  (k) , 
x  m  ’ 


k  +  k 


x 


the  finite  contribution  of  these  singularities  to  the  inversion  integrals 

constitutes  the  radiated,  freely  propagating  field.  (Singularities  occur 

also  along  the  imaginary  k  axis  associated  with  solutions  of  the  nomial- 

X  2 

mode  equation  for  negative  k  .  These  constitute  the  nonradiated, 
exponentially-damped  field  which  is  the  generalization,  for  stratified 
flow,  of  the  local  potential  flow  around  the  source.)  Using  cm  to  denote 
the  scalar  phase  speed  wm/k,  we  note  that 


k 


x 


k 


» 


which  indicates  that  no  waves  for  which  c  >  v  can  be  part  of  the  field. 

m 

It  can  be  shown  that  c^  is  always  a  decreasing  function  of  k,  so  that 
the  resonance  loci  have  two  possible  forms,  as  shown  in  Figure  3-1: 

when  v  >  c  (k  =  0)  the  source  is  "super-Froude"  with  respect  to  the  in 

m  -1 
mode  and  the  locus  leaves  the  origin  at  an  angle  of  sin  (c  (0)/v); 

*  m 

when  v  <  c  (0)  the  locus  crosses  the  k  -axis  at  a  finite  value  corre- 
m  x 

spending  to  c  (k  )  **  v,  giving  rise  to  a  transverse  wave  system.  While 
there  is  no  distinction  between  these  two  cases  in  the  formalism  to 
follow  the  inclusion  of  transverse  waves  complicates  the  numerical  pro¬ 
cedures,  and  the  numerical  algorithm  has  been  configured  for  the  super- 
Froude  case  only. 

To  evaluate  the  radiated  contributions  we  choose  the  complex  integration 
contour  of  (3-18)  to  lie  below  the  real  k  axis  so  that  the  field  will 

X 
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(CYCLES/KM) 


R-2072(U) 


u 


0 


L.J 


u 


u, 


1.1 


1  ! 


vanish  for  x  >  0,  and  we  include  the  dependence  of  u  on  k  via  k  (at 
*  r  m  x 


fixed  ky)  in  the  computation  of  the  residues, 


[k  <-2  -  4?]  _ 


dw 


=  2wv  -  2u) 


m  dk 


m  dk  dk 


to  =  ±0) 


x 


m 


=  +2(i)  v(l  -  c  c  /v  ) 
m  m  gm 


where  c  =  dm  /dk  (<c.  )  is  the  scalar  group  velocity.  The  result  is 

om  m  in 


gm  m 


5  ("^  J  ^  »  2  ) 


iV 


2v 


s  ■=>  - 


c  c  /v2)  1(J)  (z)  (f) ’  (z  )  e  +  i(%x/v 
m  gm  Ym  Ym  o  * 


(3-21) 


the  sum  carried  out  over  both  signs, 


The  numerical  procedure  is  based  on  the  observation  that  each  term 


'V, 


t,  of  the  above  sum  for  the  partial  transform  depends  on  x  only  in  its 

Itl 


complex  phase, 


v>  -  w>  e~iWmX/Vi 


(3-22) 


the  source  transform. 


(3-23) 
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can  be  computed  and  tabulated  as  a  function  of  ky  from  the  dispersion 
data  and  eigenfunction  amplitudes.  A  corresponding  table  of  complex 
phases 

-ito  Ax/v 
e  m 


can  be  similarly,  prepared,  after  which  each  £  .  can  be  generated  recur¬ 
sively  by  repeated  complex  multiplications, 


4m(x 

m 


Ax,fc  ,b) 


?  (x,k  ,z)  e 
*  y* 


-twmAx/v. 


(3-24) 


at  each  step  the  sum 

5(x,ky,z)  »  Re  |  ^  ^(x.ky.z))  ,  (3-25) 

v  m  / 

briefly  tabulated  as  a  function  of  ky,  can  be  Fourier-transformed  to 
produce  ?(x,y,z). 

Field  quantities  other  than  displacement,  such  as  scalar  strain, 
strain  rate,  and  horizontal  velocity,  can  be  deduced  from  the  simple 
relations  among  the  Fourier  amplitudes  derivable  from  the  dynamical 
equations  (3-1  through  3-4): 


i 
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(3-26) 


strain. 


strain  rate, 


■  -iuc* 


oik  _ 

velocity,  u  «  — f- 

X  kZ 


u  -  — C,' 

y  k2 


For  completeness,  the  appropriately  modified  formulas  for  the  source 
transforms  are  listed  below: 


Quantity 


m 


~  c\(l  -  c  c  /v2)  1  <f>  (z)  <J>' (z  ) 
v  m  m  gra  m  mo 


-  ~  c3k(l  -  c  c  /v2)  1  (J>^( z)  <^(zj 
v  m  m  gm  m  m  ~ 


m  m  o 


“  cfk2(l  -  c  c  /v2)"1  <P^(z)  <i>^(zj 
”  “  m  gm  m  " 


v  m 


m  mo 


(3-27) 


x 


•^r  c\(l  -  C  C  /v2)  V  (z)  <J> *  (z  ) 
m  m  gm  m  mo 


‘ik  ’ 
k 


—  C^(l  -  C  C  /v2)  (J)'  (z)  <J>'  (z  ) 

v  mv  m  gm  rm  Ym  o 


The  bracketed  quantity  [ik^/k]  is  held  outside  the  sum  in  Equation  3-25. 
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SECTION  4.  WAKE- GENERATED  WAVES 


4.1  INERT  MIXED  WAKE 


The  mixing  that  occurs  within  the  momentumless  wake  of  a  self-propelled 
body  tends  to  establish  a  more  nearly  uniform  density  in  the  wake  than  in 
the  surrounding  stratified  fluid.  Internal  waves  are  radiated  as  the 
mixed  wake  collapses  slightly  to  bring  the  internal  and  external  density 
surfaces  into  level.  For  radiated  components  whose  wavelength  is  longer 
than  the  wake  diameter,  the  effect  of  collapse  can  be  reasonably  well 
represented  by  a  time-dependent  displacement  quadrupole  in  the  crossplane. 

To  incorporate  wake  radiation  into  the  present  transform  algorithm, 
one  approach  would  be  to  model  the  wake  collapse  separately,  either 
analytically  or  numerically  [7,8]  and  use  the  corresponding  quadrupole 
profiles  as  an  additional  source  term. 

A  much  simpler,  though  not  yet  validated,  approach  has  been  used 
in  the  present  version  of  XMODE.  It  uses  a  single  parameter  in  the  form 
of  an  initial  quadrupole  strength,  and  virtually  generates  the  entire 
quadrupole  history  in  a  self-consistent  dynamical  way  as  the  radiation 
field  develops. 

The  inert  mixed  wake,  in  its  simplest  conception,  is  a  tube  of  fluid 
which,  at  a  suitable  short  distance  behind  the  self-propelled  body,  has 
an  anomalous  density  gradient  within  an  approximately  circular  cross- 
section,  no  mean  flow  perpendicular  to  the  axis,  and  in  which  turbulent 
stresses  can  be  neglected  in  comparison  to  buoyant  forces  and  inertia 
in  the  equations  of  crossnlane  flow.  This  last  assumption  means  that 
the  dynamics  of  collapse  are  adequately  described  by  the  homogeneous, 
equations  of  fluid  motion,  both  inside  and  outside  the  mixed  region. 


We  assume  that  the  only  affect  of  mixing  is  to  redistribute  density 
and  we  define  a  displacement  function  to  describe  the  redistribution, 

Ac(x,y,z-2Q) , 

such  that  if  the  mixing  were  abruptly  terminated  at  position  x,  a  fluid 

particle  on  the  streamline  originally  at  level  z  would  rise  or  fall  to  the 

new  equilibrium  level  z+A?.  Thus,  if  the  streamline  is  at  the  level  z+£, 

2 

the  particle  experiences  an  acceleration  due  to  buoyancy  -N  (C-AQ.  The 
function  A?  vanishes  for  x>o,  changes  rapidly  in  a  mixing  interval  -X<x<o, 
then  remains  independent  of  x  for  x<~X.  The  vertical  momentum  equation 
(3-1)  is  accordingly  modified  to  read 

$  +  f£  +  N2(z)  (S-A o  -  0,  (4-1) 


and  following  the  derivation  of  Section  3  with  the  volume  source  s  set  to 
zero,  we  get 


k2+j; 

\  3z" 


..  ?2  2  2 


O'  a?  » 


(4-2) 


”2  2 

where  again  the  steady  solution  obeys  £  =  v  (8/3x)  C.  For  simplicity 
we  assume  that  the  turbulent  redistribution  of  density  occurs  rapidly 
enough  compared  to  the  time  scales  x}  to  be  approximated  by  impul¬ 
sive  change, 


A? 


0,  x  >  0 


My.  z  -  zQ),  x  <  o, 


which  gives  the  source  term  in  (4-2)  the  form 


N2jj 


2 

n(x)  f  *  (x)A 

a„2 


(4-3) 


32 


l! 


i  i 

'  i 
» J 


i 


;  l 

U 


1 


i 

i 


si 

0 

n 

ii 


U 


where  n(x)  is  unity  for  x  £  0  and  zero  for  x  >  0.  The  two-dimensional 
Fourier  transform  of  this  equation  is 


or, 


+  k 


\  2 
\to 


lU.-S 

/J  ik 


V 


^(k  ,  z  -  z). 


ik  of 

x 


(4-4) 


In  terms  of  the  Green's  Function  in  (3-12),  the  solution  is 


C  = 


ik  to 
x 


2 /gm(z-  2i,i(V  zi  -  V  d2i 


(o  2 


£  -j - f  *m(z)  W* 


(4-5) 


ik  m  w  “  u 
x  m  m 


with 


w-f 


(z)N' 

m 


'  (z)A (k 


"  Zo) 


dz . 


(4-6) 


Now  for  a  range  of  wavenumbers  k  and  mode  numbers  m  such  that  the  spatial 
scales  of  the  plane-wave  normal  modes  are  much  larger  than  the  wake  radius  R, 


4f  * 
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k  R  <  1, 

y 


mR  <  Z  , 
a’ 


the  coefficients  above  are  approximately 


Qm  55  “  z)‘‘(y’z  “  Zo}  dydz 


WN  (z0^‘ 


(4-7) 


where  the  integral  Q  defines  the  wake  quadrupole  moment.  This  approxi¬ 
mation  assumes  also  that  the  ambient  Vaisala  frequency  is  constant  over 
the  wake  height  and  that  the  wake  is  axisymmetric  so  that 


If * 


dxdy  «  0. 


The  solution  transform  is  thus 


(4-8) 


The  partial  inversion  of  Section  3.3  yields  the  crosstrack  transform 


C(x,ky>z) 


|  QN2(z0)  (1  - 


/  V1 
c  c  /v  ) 
gm  m 


$m(z) 

m 


♦  '(a  ) 
m  o 


x/v 


(4-9) 
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the  sum  again  including  both  signs  of  ioj^.  This  sum,  representing  the 
raaiated  field,,  does  not  include  the  residue  of  the  singularity  at  k  «  0 

X 

which  merely  reproduces  the  impulsive  density  jump  A?  inside  the  wake  at  x 


The  source  transforms  T  generated  by  the  quadrupole  wake  are 


accordingly 


Quantity 


m 


QN^l  -  c  c  /v2)  1  A  (z)  $'(z’) 
o  gm  m  m  m 


-QN2(1  ~  c  c  /v2)  V  (z)  <J>'(z  ) 
x  o  gm  m  Ymv  Ym  o 


-QN2ia)_(l  -  c  c  /vt,)“1  $'(z)  <j> '  (z  ) 
ora  «  a  m  m  o 


(4-10) 


QN2  ~  (1  -  c  c  /v2)  ^  <p'  (z)  (z  ) 
^  o  v  gm  m  Ym  Ym  o 


-QN 


ik 

k 


ie  (1  -  c  c  /v2)  1  <p'  (z)  <p'(z  ) 
m  gia  m  m  mo 


4 . 2  SIZING  THE  QUADRUPOLE  MOMENT 

For  a  considerable  range  of  horizontal  and  vertical  wavenumber  the 
radiation  emanating  from  a  collapsing  wake  is  seen  to  be  determined  by  a 
single  constant,  the  kinematic  quadrupole  moment 


Q 


z)  A(y,  z 


zq)  dzdy. 


(4-11) 
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The  magnitude  of  this  quantity  depends  both  on  the  cross  section  of  the 
mixed  wake  and  on  the  degree  of  mixing,  as  reflected  by  the  density  level 
displacements  A.  There  is  probably  no  simple  substitute  for  a  complete 
finite-difference  simulation  of  the  turbulent  dynamics  of  wake  growth 
and  collapse  in  stratified  media  if  accurate  values  of  Q  are  essential. 
However,  certain  bounding  approximations  can  be  made  in  anticipation  of 
these  more  accurate  simulations,  and  a  quadrupole  formula  based  on  these 
approximations  has  been  included  in  the  XMODE  algorithm. 

The  estimate  derived  below  is  based  in  an  elementary  way  on  the 
dynamical  picture  of  a  growing,  momentumless  wake  given  by  Ko  [9],  with 
constants  drawn  from  the  measurements  of  Naudascher  and  Gran  [10,11]. 

The  procedure  will  be  to  assume  that  in  the  initial  growth  stage  the  density 
is  turbulently  diffused  like  a  passive  variable,  and  that  the  consequent 
growth  in  the  quadrupole  moment  continues  until  the  effects  of  stratifi¬ 
cation  intervene  rather  abruptly  at  about  one-fourth  of  the  local  Vaisala 
period.  In  Ko's  account  of  non-stratified  growth,  the  wake  radius  R 
and  turbulent  velocity  scale  u!  depend  on  time,  or  downtrack  coordinate 
x,  according  to 


R  ^  x 


1/4 


u'  'V  X 


*3/4 


(4-12) 


with  a  Reynolds  stress  that  is  given  by  an  eddy  diffusivity 

e  «  K  u'R  (4-13) 

e 

that,  is  constant  across  the  wake.  In  Gran's  measurements  this  assumption 
is  verified,  although  the  exponents  in  (4-12)  are  modified  very  slightly 
by  the  effects  of  propeller-induced  swirl.  Gran's  measurements  also 
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verify  another  assumption  of  Ko's,  that  the  wake  entrainment  rate 
R  =  vdR/dx  is  proportional  to  u', 

R  «  Kju' .  (4-14) 

To  define  the  density  distribution  inside  the  growing  wake  we  will  neglect 
the  effect  of  swirl  and  we  will  suppose,  following  a  Reynolds'  analogy, 
that  the  diffusion  of  the  passive  quantity  p  is  governed  by  the  eddy 
diffusivity  e, 


ux|^+ V‘(eVp)  =  0  (4-15) 

with  a  boundary  condition  that  equates  the  normal  flux  of  p  across  the 
expanding  wake  perimeter  at  R: 


3P, 

piR+C‘£T  =  PoR>  (4~16) 

where  u  is  the  mean  longitudinal  flow  and  the  subscripts  (i,o)  mean 

A 

inside  and  outside  the  wake  perimeter.  The  outside  density  profile  is 
assumed  linear, 


Po  “  po(zo)  ~  S(z  ~  z0)-  (4-17) 

These  equations  have  the  remarkably  simple  solution 

Pi  "  P0(z0)  “  ct(z  -  *0),  (4-18) 


as  pointed  out  by  Fernandez  [12],  since  e  depends  only  on  x  but 

3p^/9x  -  0  so  that  (4-15)  is  satisfied,  while  the  boundary  condition  becomes 


a(Rcos0)  KjU*  +  X^u'Rfacose)  =  $(Rcos0)K^u' 

where  0  is  the  polar  angle  with  z  -  zq  =  RcosQ ,  the  condition  is  everywhere 
satisfied  when 


a(Kx  +  K£)  -  BK^  (4-19) 

This  argument  suggests  that  the  growing,  non-radiating  wake  maintains  a 
constant  internal  density  gradient  whose  value  is  smaller  than  the  external 
gradient  by  an  amount  that  depends  quite  reasonably  on  the  ratio  of 
entrainment  constant  to  eddy  diffusion  constant  K£. 

The  displacements  A  in  density  surface  levels  inside  the  wake  are 
immediately  given  by 


a(z  -  zq)  =  B[(z  -  zq)  +  A], 


or 


A  = 


(z  ~  zQ). 


(4-20) 


The  quantity  (8  -  a)/8  =  E  is  a  measure  of  the  mixing  efficiency, 


(4-21) 


such  that  E  tends  to  unity  for  K  >>  K  and  to  zero  for  K,  >>  K  . 

el  i  e 
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Ti.o  kinematic  quadrupole  moment  (4-11)  is  readily  found  to  he 


Q  = 


(4-22) 


1/4 

which  is  seen  to  grow  linearly  with  time  in  view  of  the  x  -dependence 
of  R,  The  value  for  R  taken  from  Gran's  results  for  a  streamlined 
body-and-propeller  model  of  radius  yields 


R 


4 


where  RQ  =  1.19  is  the  value  at  x  =  12  R^  and  xQ  is  a  virtual  origin 
a  short  distance  (2.6  R^)  in  front  of  the  propeller.  Assuming  that 
quadrupole  growth  terminates  abruptly  after  some  fraction  f  of  the 
Vaisala  period,  when 


x  -  x 

o 


we  get  for  the  maximum  quadrupole 


(4-23) 


This  is  the  formula  implemented  in  the  SOURCE  module,  with  the  mixing 
efficiency  E  as  a  variable  input  parameter,  and  with  the  number  f  taken  as 
0.25.  The  appropriate  value  of  E  can  be  inferred  from  formula  (4-21); 
values  of  and  taken  by  Ko  from  Naudascher's  measurements,  and 
computed  by  Gran  from  his  own  data  [13],  yield  roughly  similar  values 
of  Et 


39 


Gran 

Naudascher 

Ke 

0.16 

0.18 

K1 

2.3 

2.25 

E 

0.064 

0.074 

Both  values  .of  E  are  small  enough  to  suggest  that  the  actual  mixing  tha 
occurs  in  a  propeller  wake  may  be  dominated,  by  gross  convection  due  to 
swirl.  Gran  remarks  that  a  core  of  the  wake  extending  out  to  -  0.25  R 
remains  in  solid  body  rotation  and  attains  0.8  revolutions  after  20  body 
diameters.  If  the  portion  out  to  0.5  R  were  completely  mixed,  the 

4 

equivalent  added  mixing  efficiency  would  be  on  the  order  of  (0.5)  = 

0.0625. 


4.3  TREATMENT  OF  HIGHER  MODES 


The  approximation  (4-7)  for  the  mode-dependent  quadrupole 
coefficients  (4-6)  overemphasizes  tha  higher  modes  and  higher  wave- 
numbers  by  treating  the  source  like  a  point  quadrupole  singularly  in 
the  crossplane.  The  normalized  mode  product  appearing  in  the  source 

transform  formulas  for  strain,  velocity,  and  strain  rate  is  on  the 
average  an  increasing  function  of  m,  of  order 


_2  f, 2  ,  Pm 

$  4  ^  P  N>  dz  =  —z 
YmYm  2 

where  one  defines  a  vertical  wavenumber  p  by 

rm 

dz  =  dz 

■  m  J  Tm  J  Ym 


ft' 

Jr  m 


and  a  mode-averaged  N  by 


H2  /V  dz  =  A2  N2dz  »  1. 
my  m  J  m 
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One  can  easily  derive  from  the  normal-mode  equation  (3-10)  the  relation 


.  2  .  2  -2  2 
K  +  p  =  c  N  , 
rm  mm 


(4-24) 


so  that  the  wake-generated  source  transforms  for  strain, 

T(wake)  *  QlTc"  (1  - 

Nm 


tend  to  grow  without  mode  limit.  The  extra  factor  c^k  in  the  body¬ 
generated  source  transforms  attenuates  these  functions  at  higher  mode 
numbers.  This  behavior  is  evident  in  the  body-dominated  and  wake-dominated 
cases  illustrated  in  Figures  2-5  and  2-6  . 

The  actual  wake  excitation  coefficients  defined  in  equation  (4-6) 
will  diminish  with  increasing  k  and  m  at  a  point  where  the  scale  of  the 
mode  function  becomes  comparable  to  the  maximum  wake  radius  R.  A  simple 
form  factor  to  simulate  this  effect  qualitatively  has  been  included  in 
the  wake  transform  algorithm, 


2  „2  ’ 


1  +  0.4  q  R' 


where  q  is  a  total  wavenumber  defined  by 


2  ,2,2 
q  =  k  +  P  , 
’m  m 


(4-25) 


(4-26) 


and  computed  very  simply  from  quantities  on  hand  by  formula  (4-24)  above. 
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